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Abstract 

>'. 

q , We construct a local realistic hidden-variable model that describes the states 

^ | | and dynamics of bulk-ensemble NMR information processing up to about 12 

qq ■ nuclear spins. The existence of such a model rules out violation of any Bell 

CN ■ inequality, temporal or otherwise, in present high-temperature, liquid-state 

. NMR experiments. The model does not provide an efficient description in 

that the number of hidden variables grows exponentially with the number of 
^ | nuclear spins. 
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High-temperature, liquid-state nuclear magnetic resonance (NMR) provides a testing 
q ■ ground for the new ideas for information processing that are being developed in quantum 
information science The fundamental information-processing elements used in NMR 
Qh! are two-level nuclear spins, called qubits, which are bound together in a single molecule. A 
liquid NMR sample contains a macroscopic number of molecules, each of which functions 
as an independent information-processing unit. The molecules are initially in thermal equi- 
librium at high enough temperature that the nuclear spins are only weakly polarized along 
the direction of a strong magnetic field. NMR techniques cannot control the quantum states 
of individual molecules, and the measurements performed in NMR detect the average mag- 
netization of the entire sample. For these reasons the use of high-temperature, liquid-state 
NMR to emulate quantum computation is called bulk- ensemble quantum computation. 

The original proposals for quantum information processing using NMR were greeted 
with enthusiasm tempered by skepticism. The enthusiasm led to a remarkable series of exper- 
iments in which NMR techniques have been used to implement the operations for a variety 
of quantum-information-processing jobs involving up to seven qubits (for reviews of NMR 
information processing, see Refs. The persistent skepticism has to do with questions 

about the "quantumness" of NMR information processing. Initially based on doubt that the 
highly mixed states used in NMR could be used to achieve genuinely quantum-mechanical 
effects, these questions were made concrete by the realization that all the quantum states 
accessed in present experiments are unentangled ||10|| . Entanglement is often thought to be 
an essential feature of quantum computation pd|,|T2f. Arguments for an essential "quantum- 



ness" in NMR information processing are presented in Ref. ||13|| , and an entirely different 
method for characterizing the "quantumness" of NMR is developed by Poulin [[HJ . 

The absence of entanglement in present NMR experiments means that the statistics of 
measurements made at any time during the experiments can be understood in terms of a 
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local realistic hidden-variable model |15| in which each spin has objective properties that de- 
termine the results of the measurements. In a local realistic model, the correlations observed 
in experiments can be attributed to classical correlations between realistic properties of the 
component qubits. The ability to describe the correlations observed in NMR experiments 
in terms of classical correlations between the qubits casts doubt on the "quantumness" of 
the experiments. 

Modelling the statistics of the measurements made in NMR experiments is not sufficient 
for understanding the experiments. One must also be able to model the dynamics of the 
nuclear spins. NMR experimenters can implement with high accuracy any unitary opera- 
tion, including the nonfactorizable unitary operations — those that cannot be written as a 
product of unitaries for each qubit — that produce entanglement when applied to pure quan- 
tum states. Previous attempts [16j to devise a local realistic description of the dynamics 



were only partially successful in that they did not provide a local realistic description of 
the changes produced by nonfactorizable unitaries which reproduced all the predictions of 
quantum mechanics. This left open the possibility that one might not be able to describe 
the correlations observed in successive measurements separated by nonfactorizable unitary 
operations in terms of local realistic properties and thus that present NMR experiments 
might violate temporal Bell inequalities W7\ for successive measurements. 



In this Letter we report a local realistic hidden-variable (LRHV) model for the states 
and dynamics of bulk-ensemble NMR information processing up to about 12 qubits. The 
existence of such a model rules out violation of any Bell- type inequality in present NMR 
experiments. This conclusion applies only to the bulk-ensemble model of information pro- 
cessing realized in present high-temperature, liquid-state NMR experiments; it does not 
apply to NMR methods based on distilling a pure state from a thermal state [[5,18,19 . Our 



model is not satisfactory from the point of view developed by Schack and Caves jTjJ because 
the number of hidden variables scales exponentially with the number of qubits. 

All NMR quantum computing experiments performed so far work in the following way. 
The state of each molecule, consisting of A" active spin-| nuclei, is described by a density 
operator 

p=(l-e)i/2 N + ep 1: (1) 

which is a mixture of the desired state of the quantum computer, pi, with the maximally 
mixed state for A" qubits, 1/2 N , 1 being the unit operator. When p\ is a pure state, p is 
called a pseudopure state 0. 

The molecules in an NMR sample begin in thermal equilibrium, with a weak polarization 
a = hv/2kT ~ 2 x 10~ 5 at room temperature, where v ~ 300 MHz is the average resonant 
frequency of the active spins in the strong longitudinal magnetic field. The first step in NMR 
information processing is to transform the molecules from equilibrium to a pseudopure state 
PJ2H|] . A consequence of pseudopure state synthesis is that the mixing parameter scales like 
e = aN/2 N . 

After synthesis of the desired initial state, the computation begins. The unitary oper- 
ations required for the computation can be constructed from sequences of radio-frequency 
pulses alternating with periods of continuous evolution under the nuclear-spin Hamiltonian 
A unitary operator U takes an input state p to an output state 

Up l/t = (1 - e)l/2 N + eU Pl U^ . (2) 
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The maximally mixed state is unaffected by the unitary transformation. The output state 
retains the form ([]]) with the same value of e, and — this is the essence of the bulk-ensemble 
paradigm for quantum computation — pi undergoes the desired unitary transformation. 

The computation completed, the last step is to read out the answer. By applying radio- 
frequency pulses and then measuring the transverse magnetization of the sample, an NMR 
experimenter can determine the expectation value of any product of spin components, one 
for each qubit 0-0. These expectation values have the form 

C(a) = C(a 1 ,...,a N ) 

= tr(p <r • ai ® • • • ® <r • ajv) 

= etr(pi <x • ai Cg) • • - ® & ■ ajv) • (3) 

Here and throughout, a tilde over a quantity denotes a collection of N such quantities, one 
for each spin. In Eq. (^) the tensor product includes one operator for each spin; the vector 
operator cr = le + <J x e x + <J y e. y + a z e z , where a x , a y , and a z are the Pauli operators; and 
a r is either a spatial unit vector, in which case er • a r is the component of spin r along the 
direction a r , or the unit vector eo in the "zero" direction, in which case spin r does not 
contribute to the expectation value. The last equality in Eq. ([|) assumes at least one of the 
vectors a r is a spatial direction. 

The expectation values (|3|) express the correlations between spin components of different 
spins. The maximally mixed state does not contribute to the correlation coefficients, which 
are determined by the state pi that undergoes the desired evolution. The mixing parameter 
e measures the strength of the magnetization signal. The scaling e = aN/2 N that comes 
from pseudopure state synthesis thus leads to an in-principle demand for an exponentially 
increasing number of molecules as the number of qubits increases ||2,21|, implying that bulk- 



ensemble quantum computation is not suitable for large-scale quantum computation. 

Before constructing our LRHV model, we recall that any iV-qubit density operator f has 
an associated quasidistribution 



w?(n) =tr(tQ(n)) , (4) 
where the vectors in the set n = (toi, . . . , tijv) are spatial unit vectors and 

Q(n) = -^(i + 3ni • &) <g> • • • <g> (i + 3n N ■ a) . (5) 
For each spin, the unit vector n can point in M different directions satisfying 

= 5>*, (6) 

n 

\ 5 ik = JfY, n i n k > ( 7 ) 



n 



where the sums are over the possible directions and the subscripts indicate spatial com- 



ponents of n. Condition ([?]) means that the vectors y3/jVn form a resolution of the 
3-dimensional unit tensor; condition (|6]) places an additional constraint on the placement of 
the vectors. The vertices of a tetrahedron give the minimum number, M = 4, of possible 
directions. The six vectors along the cardinal directions make up another simple possibility. 
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The density operator is then given by [p2]] 



T = Y^wr(ri)\n){n\ , (8) 

n 

where \n)(n\ = |n 1 )(n 1 | eg) • • • ® |njv)(njv| and |n)(n| = |(1 + <x ■ n) is the +1 eigenstate 
of a • n. In terms of Wf (n), the correlation coefficients (Q) take the form 

= H (n)(ai • mi) ■ • ■ (a/v ■ m N ) , (9) 

IT. 

where m r = n r + e . 

Under a unitary transformation the quasidistribution evolves according to 

w ^(h') = J2TLw,(h). (10) 



n 



Here the transformation matrix T u has matrix elements 

TL = (n\tfQ(h f )U\h) = v, mw fr(h') , (11) 

given by the quasidistribution for U\h). 

A separable density operator is one that has an ensemble decomposition in terms of prod- 
uct states. Such a state has no entanglement. If the quasidistribution wp(n) is everywhere 
nonnegative, then p is definitely separable, and the statistics of all measurements can be 
understood in terms of classical tops whose probability to point in the directions fi is w p (h). 

For any density operator, the quasidistribution satisfies Wp(h) > [minimum eigenvalue of 
Q(n)} = {-2)A N - 1 /M N = -2 2N ~ 1 /Af N pH^. Thus for density operators of the form (|), 
the quasidistribution is everywhere nonnegative if ]TU| , p2| 



^r^ s "- < 12) 

Such states are unentangleable by any unitary transformation. For the polarization a ~ 
2 x 10~ 5 of present NMR experiments, all states up to about 12 qubits are unentangleable. 
It is known P3| that entangled states of the form ([I]) exist for e > (1 + 2 JV ~ 1 )~ 1 = rj', i.e., 
N > 2/of, but whether there are entangled states for 7] < e < rj' is an open question. 

We turn now to constructing a LRHV model for unentangleable states, i.e., for e < rj. 
A straightforward model regards the directions n as hidden spin directions that determine 
the results of measurements stochastically. The problem with this simple model is that 
the change in the probability distribution w p (n) as a consequence of a unitary transforma- 
tion should be described by transition probabilities, which give the probability to go from 
initial directions n to final directions h'. The matrix elements (|TT| ) provide candidate tran- 
sition probabilities, but they cannot be interpreted as transition probabilities because they 
generally take on negative values. Schack and Caves ]I6| derived nonnegative transition 
probabilities from these matrix elements, but the dynamics described by these transition 
probabilities did not duplicate the predictions of quantum mechanics. 

To include the unitary dynamics in the local realistic description, we construct a deter- 
ministic LRHV model that includes the probabilities w p (n) in the set of hidden variables. 
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A quasidistribution Wf(h) can be regarded as a vector with Af N components labelled by 
the directions n. An arbitrary vector in this space, denoted by w, has components w(n). 
We are only interested in normalized vectors, i.e., J^nWi™) — 1- 

The hidden variables in our model are a vector w, a set of spin directions n, and a set of 
real variables A = (Ai, . . . , Ajv) such that — 1 < A r < 1, r = 1, . . . , N. The hidden variables 
are denoted collectively by A = (w, h, A). The probability density for A is 

P(\) = ±5{w-iDfi)w(n), (13) 

indicating that the variables A are distributed randomly, the hidden vector w has a definite 
value given by wp, and the probabilities for the hidden spin directions are determined by 
the hidden vector w and, hence, are written as w(h). One can think of the hidden vector 
w as a set of parameters that weight a "roulette wheel," so that when the wheel is "spun" 
to generate hidden spin directions, the probability for directions n is w(h). Notice that the 
model requires that all components of wp be nonnegative, since these components become 
probabilities for the spin directions. 

Measurement results are governed by functions A r (a r , A), r = 1, . . . , N. The value of 
A r (a r ,X), either ±1, determines the result of a projective measurement of the component 
of spin r along spatial unit vector a r . The model is realistic because the results of spin- 
component measurements are determined by the hidden variables, and it is local because 
the result of a measurement of the component of spin r along direction a r depends only on 
a r and the hidden variables, not on measurement directions for other spins. 

We choose the spin-component functions to be 

M°^)=M°r,Kn r )={ + _i < i4 > 

By using m r = n r + eo, we can employ this function in cases where spin r is not involved in 
a measurement, i.e., when a r = e , which gives A r (e , A) = 1. The correlation coefficients 
predicted by the LRHV model, 

^LRHv(a) = J d\P(X)A 1 (a 1 ,X)---A N (a N ,X) 
= ^2 J dw 5(w — Wp)w{h) 

n 

N 1 f+l 

x TT - / dA r A r (a r ,A r ,n r ) 
r=i 2 J - 1 

= J2 wp(n)(a\ ■ mi) • • ■ (a N ■ m N ) , (15) 

n 

duplicate the predictions of quantum mechanics. 

Unitary dynamics fits easily into this model. For any unitary transformation U, the 
hidden variables for each molecule are updated in the following way: the variables A are 
chosen randomly, the hidden vector w is updated deterministically using the transition 
matrix as in Eq. (|10D , and a new set of hidden spin directions is obtained by spinning a 
roulette wheel weighted by the new hidden vector, w^^. By treating the components of 
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the hidden vector not as probabilities, but rather as parameters that are used to generate 
hidden spin directions stochastically, we avoid the need for a nonnegative transition matrix. 

This description of the updating applies to a discrete unitary transformation, but it is 
easy to generalize it to a quasicontinuous hidden-variable dynamics: each molecule updates 
randomly in the fashion just described; i.e., each has a probability jdt to update within each 
time interval dt. The only requirement on this quasicontinuous dynamics is that the mean 
time 7 _1 between updates be large compared to the precession time v~ x of the nuclear spins 
in the strong magnetic field, but shorter than the duration of the radio-frequency pulses 
that are used to produce the desired dynamics. 

It is trivial to generalize the LRHV model presented here to nonunitary evolutions, since 
these evolutions are, like unitary transformations, linear in the density operator. It is likely 
that the model could be extended to include all dynamics that accesses only separable 
states, i.e., states that have an expansion like Eq. (H), but with more general nonnegative 
quasidistributions than the canonical form ([|). 

The LRHV model developed here achieves our purpose of determining whether present 
NMR experiments can violate Bell inequalities. Bell inequalities are founded on two as- 
sumptions: the assumption that systems have objective properties and a "no-disturbance" 
assumption that asserts that the relevant measurements report faithfully the values of 
these properties. Standard Bell inequalities [15] justify the no-disturbance assumption from 
locality — measurements here cannot affect properties there. Temporal Bell inequalities [jT7 



which involve successive measurements on a system, have a tougher time justifying the no- 
disturbance condition. Our LRHV model for NMR experiments describes the statistics of all 
measurements in terms of evolving classical correlations between realistic properties of the 
constituent nuclear spins. The conclusion is that NMR experiments up to about 12 qubits 
cannot violate any Bell inequality, temporal or otherwise. 

Our purpose achieved, we acknowledge that our LRHV model is terribly contrived. It 
succeeds in giving a local realistic description of the dynamics by the brute force device 
of including an encoding of the entire density operator among the hidden variables. As a 
result, it requires an exponentially increasing number of hidden variables, ~ 4^ in the most 
efficient version of the model. It leaves open the possibility that [|16[] the "quantumness" of 
NMR information processing lies in the ability to implement nonfactorizable unitaries that 
do not have an efficient local realistic description. 

This work was partly supported by the National Security Agency (NSA) and the Ad- 
vanced Research and Development Activity (ARDA) under Army Research Office (ARO) 
Contract No. DAAD19-01- 1-0648. CMC received support from National Science Founda- 
tion Contract No. PHY99-07949 at the Institute for Theoretical Physics of the University 
of California, Santa Barbara. 
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